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Under certain conditions, tiie quantum ^-kicked iiarmonic oscillator displays quantum resonances. We con- 
sider an atom-optical realization of the (5-kicked harmonic oscillator, and present a theoretical discussion of the 
quantum resonances that could be observed in such a system. Having outlined our model of the physical system 
we derive the values at which quantum resonances occur and relate these to potential experimental parameters. 
We discuss the observable effects of the quantum resonances, using the results of numerical simulations. We 
develop a physical explanation for the quantum resonances based on symmetries shared between the classi- 
cal phase space and the quantum-mechanical time evolution operator We explore the evolution of coherent 
states in the system by reformulating the dynamics in terms of a mapping over an infinite, two-dimensional set 
of coefficients, from which we derive an analytic expression for the evolution of a coherent state at quantum 
resonance. 
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I. INTRODUCTION 

As the field of quantum chaos (V. '2] has emerged over the 
last few decades, much research has been directed towards the 
quantal counterparts of classically chaotic Hamiltonian sys- 
tems |3, 4|. Amongst these systems, one-dimensional, peri- 
odically kicked quantum systems have proved a popular tar- 
get for study. The dynamics of these systems, which are de- 
scribed by discrete time evolution operators acting over the 
period of one or more kicks, often display quantum reso- 
nances and antiresonances; changes in the degree of local- 
ization of the dynamics which are strongly dependent on the 
value of an effective Planck constant. Advances in the cool- 
ing and trapping of atoms and ions has paved the way for 
substantial direct experimental investigation of such archety- 
pal quantum-chaotic systems, particularly the 5-kicked ro- 
tor [5i\g U] M M M MM M M and also the 6- 
kicked accelerator |[l6l|T7l[l8l[l9l|20l|2l]l22llll- An- 
other possibility is the ^-kicked harmonic oscillator (KHO) 
El 125] |26l 1221 |28l im |30l; many kicked atom (ion) config- 
urations involve trapping potentials which are to some degree 
harmonic, making the dynamics of the KHO a relevant con- 
sideration for a wide variety of experiments. 

The classical KHO was originally developed as a model 
for the interaction of charged particles with electromagnetic 
fields ll3n [32l [33l . The frequency degeneracy of orbits in the 
harmonic oscillator renders the KAM theorem inapplicable 
to the KHO (i.e., all invariant tori of the harmonic oscilla- 
tor may be destroyed at any finite perturbation strength) li34J . 
In cases where the kicking period rationally divides the nat- 
ural oscillator period this allows the KHO to display exotic 
non-KAM chaos, manifest as interlinked regions of chaotic 
dynamics forming stochastic web structures throughout phase 
space 1331 l35l . The quantum KHO has also been studied in 
considerable depth |28 , 3dv3X,^ 39J; particularly in the case 
where it is kicked four times per oscillator period, when it can 
be related to the kicked Harper model Il40ll4ni42ll43ll44ll45l. 
Depending on the relative alignment of the kicking potential 
with respect to the trapping potential, the KHO can display ei- 
ther quantum resonances (even kicking potentials) or quantum 



antiresonances (odd kicking potentials) Il38l . 

In this paper, we consider the effects of quantum resonances 
that may be observed in an atom-optical KHO. In Section [ll| 
we outline our model of the atom-optical KHO and introduce 



the kick-to-kick Floquet operator. In Section III we derive the 



structure of the quantum resonances and their effects on the 
dynamics. We consider the relation between the quantum res- 
onances, the eigenstates of the Floquet operator for q kicks, 
and that Floquet operator's expansion in terms of phase space 
displacement operators; this leads to an alternative, physi- 
cally intuitive explanation for the occurrence of quantum reso- 
nances. In Section|lV]we consider the measurable manifesta- 
tion of quantum resonant dynamics in the atom-optical KHO, 
including the experimental accessibility of the necessary pa- 
rameter regimes. Section[V]consists of the conclusions, which 
are then followed by three technical appendices. 



II. ATOM-OPTICAL 5-KICKED HARMONIC 
OSCILLATOR 

A. System Hamiltonian 

We consider a two-level atom (or ensemble of noninter- 
acting atoms) or ion, subject to tight radial confinement and 
relatively loose harmonic axial confinement, and periodically 
driven by pulses from an off-resonant laser standing wave. We 
assume the pulse duration to be sufficiently short for the par- 
ticle centre of mass displacement during the pulse to be negli- 
gible relative to the period of the standing wave (Raman-Nath 
regime). To a good approximation, the excited internal state 
can be eliminated, and the system can be described by the 
one-dimensional single particle Hamiltonian 



H 



32 Ma/x^ n^^tn 



2M 



8A 



cos{Kx) 6(t-nT), (1) 



describing the centre of mass motion in the axial direction 
only. Here M is the particle mass, and o) the frequency of 
the harmonic trapping potential; the kicks are pulses of dura- 
tion tp, repeated with period T, from the laser standing wave. 
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detuned by A from the two-level transition and with Rabi fre- 
quency Q, where K - AnjAi and Ai is the laser wavelength. ' 

We have set the origin of the x-coordinate to be at the mini- 
mum of the harmonic potential, and choose the laser standing 
wave to be aligned symmetrically about this minimum. The 
resulting cosine potential is an even function in Hamiltonian 
([T]l; hence, we can expect quantum resonances, i.e., particu- 
lar parameter regimes associated with ballistic growth in the 
system energy, to be observable Il28l [38l . Shifting the stand- 
ing wave by /Il/8 yields a sine potential — an odd function, 
producing antiresonances in parameter regimes where quan- 
tum resonances occur for a cosine potential. When the align- 
ment of the laser standing wave with respect to the trapping 
minimum lies between these extremes, intermediate quantum 
resonant behaviour can occur 1281 |3F |. It is therefore not 
essential for the laser standing wave to be aligned perfectly 
symmetrically in order for quantum resonant behaviour to be 
observable, although from now on we assume this, for sim- 
plicity. More important is that any kick-to-kick variations in 
the standing wave alignment be negligible compared to the 
standing wave period, otherwise the resulting randomness in- 
troduced into the dynamics will destroy the desired quantum 
resonant behaviour. 



B. Time evolution 

The time evolution of the quantum KHO is described by a 
Floquet (discrete time evolution) operator F which maps the 
state of the system at a time immediately prior to one kick 
onto the state at a time immediately prior to the next. It is 
convenient to describe the Floquet operator F in terms of the 
harmonic oscillator creation and annihilation operators a ' , a, 
where a - yfma)j2h(x + ip/mcj). Hence, 

p — g-/(fl*a+l/2)Tg-/*rcos[/;(fl+a*)]/ V2/;2 ^2) 

where r = lliT is the dimensionless kicking period, k - 
's/lAMo) is the dimensionless kicking strength, 
and - K y/h/lMo) is the Lamb-Dicke parameter, which 
characterizes the scale of the harmonic ground state relative 
to the laser wavelength. The parameters t and k emerge natu- 
rally from a rescaling to dimensionless dynamical variables of 
the KHO classical map |3T1[32][33], whereas jf now assumes 
the role of an effective Planck constant ||29| . 

C. Pliase space symmetries 

We confine our attention to when the kicking period ratio- 
nally divides the natural period of the oscillator, i.e., when 
T - 2nr/q for coprime natural numbers r, q. Except in 
the special cases where q = 1 or 2, the phase space of the 
classical KHO is then spanned by a stochastic web ||33l [35l , 



The quantity hK is tlierefore equal to two photon recoils. 



i.e., interconnected regions of chaotic dynamics penetrating 
through all of phase space. If ^ e q^ = {1,2,3,4,6) the clas- 
sical phase space has a discrete translational (crystal) symme- 
try; otherwise the stochastic web has a quasicrystal structure. 
These crystal symmetries are reflected in the quantum regime 
through the existence of infinite sets of unitary displacement 
operators D(a) — e™ " commuting with F'', the evolution 
operator for one natural oscillator period 125,26,29,^1. This 
implies that eigenstates of F'' for q e qc have a comparable 
crystal structure, and, importantly, are consequently extended 
in phase space. In the non-trivial crystal cases q € {3,4,6) 
— when the displacement operators form a two-dimensional 
crystal lattice in phase space — the set of displacements com- 
muting with F'' is {D{yq)], where 

74 - kn/T] + Htt/i], (3) 
73 = 76 = ^ (l + // V3) ;r/;7 + Z (l - // Vs) n/rj, (4) 

for arbitrary integer k, I? 

III. ORIGIN OF QUANTUM RESONANCES 
A. Quantum resonant values of ?]- 

The general study of quantum resonances in systems with 
periodic classical phase spaces has been carried out in consid- 
erable theoretical detail by Dana |28 36 37 1. Our intention 
here is to provide an accessible derivation of the conditions 
under which quantum resonances occur specifically for the 
atom-optical KHO, as described in Sec. |ll] 

With this in mind, we substitute r = 2nr/q into (|2]), and 
rearrange F^' into the product 

9-1 

pq — g!(a*S+l/2)2n-;V/?g-kcos[i7(a+flt)]/\^;;-g-i(fl*S+l/2)2n-j>/? 

(5) 

We use the identity e"'^"'^(a + a^)e - ae " + a^e" to trans- 
form each term of (jsj, and write F'' as 

9-1 

F" = Y[ (-'-^cos [77 [ae-'^^i'^i + a'^e'^"^''^")] / ^/It]^) , 

j=0 

(6) 

that is, as a sequence of q cosine kicks directed along q axes 
equally distributed in angle about the origin in phase space. 
Using the identity [47| 

^*cos© ^ ^ ijjj(;^)e-j^ (7) 
j=-°° 



^ Here a is the conventional complex index of a coherent state, such that 

D(a)\0) = \a) anda|a> = a\o). 
^ To reinforce the connection to the classical stochastic web, one can rescale 

the classical realization in terms of an appropriate dimensionless position 

- momentum pair and obtain (for the cases q = 3,4, 6) the translations J3} 

and |4| with set to unity [30j. 



3 



(the Jj are y'th-order cylindrical Bessel functions), we now ex- 
pand F'' in terms of products of unitary displacement opera- 
tors: 

q-\ DO 
;=0 kj=-oo 

Taking the commutation relation for the unitary displacement 
operators D(a)D(fil) = e"P' 1^ b(Ji)b{a) El, it follows that 
for a = ik^jje'^"""'"', /? = ik„7]e'^'""''i , then e"'^'^"''^ = 1 (i.e., 
the displacement operators commute) whenever 



T 2, , • l2nr[m-n]\ 
2n k,„k„ sm = Zttw, 

q 



(9) 



for integer w. If (j9| is fulfilled for all m - n e {0, 1, . . . , <7 - 1) 
it follows that all displacement operators in ( 8| mutually com- 
mute, and hence all of the product terms in (6]l must mutually 
commute. Hence, under these conditions, the vth power of F'^ 
will describe a ^-fold diffraction exactly equivalent to (|6]l but 
with amplified kicking strength kv: 

9-1 

j=o 

(10) 

This is analogous to the case where the quantum resonance 
condition is fulfilled for the quantum (5-kicked rotor, when the 
free evolution operator between kicks collapses to the iden- 
tity, and the evolution induced by v iterations of the Floquet 
operator is equivalent to that caused by a single kick with v 
times the kicking strength f/l ^ '46| . 

The values of jf for which (|9j is fulfilled therefore deter- 
mine when quantum resonances occur in the KHO. In sum- 
mary, (jojl is fulfilled for: all values of 77^ when q e {1,2); 77^ 
equal to a multiple of n when q - 4; and rf equal to a mul- 
tiple of In/ V3 when q € {3, 6) (see Appendix |a|i. For q t qc 
it is impossible to satisfy (|9]l. Note that, for ^ e {1, 2), (j6]l re- 
duces to being exactly equivalent to what one would expect for 
the quantum 5-kicked rotor fulfilling the conditions necessary 
for quantum resonance; the free evolution between kicks has 
no net effect, and the (time-periodic) dynamics can be equiv- 
alently described by a single amplified kick. However, the 
fact that the free evolution between kicks has no net efifect is 
equally true classically (for q = 1 this simply means that after 
one oscillator period between kicks, the oscillator has returned 
to its original location in phase space), and one must therefore 
conclude that quantum resonances do not occur for q e {1,2), 
as is consistent with the dynamics' independence of the effec- 
tive Planck constant 77^. 

Consequently the KHO displays quantum resonances only 
if ^ 6 {3,4,6), and for 77^ equal to multiples of the princi- 
pal quantum resonance values: n when q - 4, and 2n/ V3 
when q - 3 or q = 6. At further rational multiples a/b of the 
principal quantum resonance values, commutativity between 
kicks depends on the value of k,„k„ inj9]), and only occurs for 
some terms within the summation in ([8]), generating higher or- 
der quantum resonances. The proportion of terms for which 
the displacements commute decreases with the denominator 
b, approaching zero as a/b approaches irrationality. 



B. Relation to phase space symmetries of 

When q e {3,4,6), the set of unitary displacements in ([8]) 
describe a two-dime nsion al crystal lattice in phase space. In 



analogy with Section II C we denote this set {DiFq)], where 



r4 = A;77 -H 7777, 



rv3 



+ 1 



77 + 1 



rv3- 



(11) 

(12) 



Comparing ( fTT) and ( fTS] ) with Q and Q, one sees that, at 
the principal quantum resonance values [77^ — n (q — 4), or 
77^ = 27r/ V3 (q - 3,6)], the set {D(F^)) of displacements 
constituting F'' is identical to the set {Diyq)} of translational 
symmetries of the eigenstates of F^. Moreover, for 77^ equal to 
rational multiples a/b of these principal values, D(Tq{k, I)) re- 
mains an element of [Dijc^] whenever k and I are integer mul- 
tiples of b. This (complete or partial) equivalence of {D{Yq)] 
and [Diyq)] offers a physical explanation for quantum reso- 
nances in the KHO; at quantum resonance values of 77^ the 
displacements constituting each kick are translational symme- 
tries of F^'s eigenstates, and consequently commute with F^. 



C. Evolution of an initial coherent state 

I. Motivation 

In the context of cold, trapped atoms or ions, the initial mo- 
tional state of the atom(s) or ion in any experiment will gener- 
ally be as close to a zero-temperature state as is practicable for 
the specific experimental configuration. This provides a clear 
motivation to study the system evolution when the initial con- 
dition is a harmonic oscillator ground state. In the case of an 
initial atomic Bose-Einstein condensate, this would describe 
the initial condensate mode if atom-atom interactions could 
be neglected, which could potentially be achieved with the 
aid of magnetic fields to exploit a Feshbach resonance |49|. 

Coherent states are harmonic oscillator ground states which 
have been displaced in phase space, and as such can in prin- 
ciple be accessed by displacing an initial ground state in mo- 
mentum (e.g., by application of a Bragg pulse) or position 
(e.g., by suddenly changing the trap minimum), possibly in 
combination with a free evolution subsequent to the displace- 
ment. In what follows we therefore consider the evolution 
of an initial coherent state analytically, which, in the case of 
quantum resonance, proves to be comparatively tractable. 



2. Phase space lattice state expansion 

We begin by applying the identity (jTjl to the Floquet opera- 
tor yielding 

00 

P ^ ^-i{a^MI2)2nrlq J ^{^)D (ir^k) , (13) 



4 



where we have introduced ^ = -kI V277^ = —Q^tp/8A It 
follows that applying F to an arbitrary initial coherent state 
\a} will yield a weighted superposition of regularly spaced co- 
herent states. More generally, we will show that if the state of 
the system \iffj} after j applications of ( pjj i takes the form of a 
"lattice state" 

oo 

m,n--oo 

(14) 

then \t//j+\} = f'l'/')) will take an essentially similar form. The 
\ae are all coherent states, and the M|,f^„ a set of 



complex coefficients 



Applying (13i to (14i, and using the identities 



b{a)I)iJ5) = e^''P'-"'P^l^b(a+li) and e-'"' t) (a)\l3) = 
D{ae^''^)\Pe^''^), now yields, in the first instance. 



k,m,rf=-oo 

X b {iri [(m + k)e-'^'''li + 

X e"'(-'+l)''/9|Q,g-'2TO'+l)''/9^ 



(15) 



When q e {3, 4, 6), e^^'^'"'!'' can be written as an integer super- 
position of unity and g^'^'"'/?^ allowing ( 15 1 to be simplified. 
From now on we assume r = 1 to illustrate the procedure: 



_l_e-'-2;r/9 if^ = 3, 
-1 if? = 4, 

-l+e-'^^li if? = 6. 



(16) 



Using ( 16 1 to make replacements in ( 15 1, one obtains 



k,m,n--co 

X b {iri [-n + (m + n^q + k)e'''^''^''^ 
where we have introduced 



(17) 



-1 if? = 3, 

if? = 4, 

1 if? = 6. 



(18) 



Effecting a change of indices such that m ^ m' - -n and 
n ^ n' - m + n^g + k now gives 



2 i'j^iOM]'' 



^-ikmrj- sin(2;r/(j) 



k,m,n=-oo 

X b [ii] [m + ne-''2'^/«])e-'<-'+i>/'?|ae-'2''(-'+i)/?)^ 



(19) 



where we have dropped the primes, for convenience. With 



(19 1 we have an expression of the same form as ( 14 1, but for 



f will usually be a positive quantity, as this corresponds to red detuning. 



the state and can obtain the evolution resulting from F 

as a mapping of the lattice coefficients: 



[j] ^-iktmf sin{2n/q) 

m^^,+n-k,-m 



(20) 



The coefficient mapping ( 20 1 fully describes the time evo- 
lution of lattice states in general. The coherent state at the 
center of the lattice state when j - mod ? can be chosen 
arbitrarily. Hence, choosing an initial set of coefficients 



(21) 



amounts to following the evolution of an initial coherent state 
\a). 

The general evolution of the coefficients over many kicks 
described by ( [20| yields an ever-lengthening product of nested 
sums that is not especially insightful. In the case of quantum 
resonance, however, t;^ sin(27r/?) - wn, which simplifies the 



final exponent in ( 20 1, allowing a simpler expression for evo- 
lution over many kicks to be extracted. In the next section we 
follow the evolution of an initial coherent state, for ? = 4 and 



77 = n, by repeatedly applying the mapping ( 20 1 



3. Analytic solution at quantum resonance when q = 4 



Our initial condition is a single coherent state, as described 
by (21 1. We first rearrange the indices to obtain m'j°\_,„ = 
Sn-kfiO-„,fi = S„^ii6„,fl, which we substitute into (20 1, yielding 



Mill = Z i'm)Sn.k6„,,oe-"""'' = ("UOSm.oe-'"'"". (22) 



We now repeat this procedure; rearranging the indices and us- 
ing the identity i"J„(0 - i "J-n(0 yields 



n—k,-m 

which we substitute into (|20li: 



M 



[2] 

n-k,-m 



i{n-k)mn 



(23) 



(24) 
(25) 



We next substitute ( [25| ) into ( [20) 1, accompanied by a change of 
indices k ^ -k and use of i''Jk{0 - i '^J-kiO'- 

00 

<l = f'JJire"""" JkiOJn.kiOe"'". (26) 

/;=-oo 

The summation over k can be reduced using Graf's addition 
theorem (Appendix [B| to give 



(27) 



We repeat the process, again with a change of indices k ^ -k 
followed by the application of Graf's theorem: 

Ml^k.-,n = i"-'jn-k(0i"'J,n(20e-'^"-''^'"'', 



(28) 



k=-oo 
imnn 



(29) 
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FIG. 1: (color online). Phases of the coefficients m]!^} as given by 
^30\ , which describe the state of a KHO (with q = 4) after the 
application of N (N > 2) kicks to an initial coherent state, for the 
quantum resonant case rf = n. The phases shown neglect the par- 
ity of the Bessel functions in the coefficient [that is, we actually plot 
M'v! /■/,., (C|"V)-/«(Cj"Y) = (-l)™!""^"] so that the pattern remains 
constant for all N > 2. 
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FIG. 2: (color online). Q functions Qj(a,a') = j\a)\^ I n f48l for 
the KHO (with q = A and k = -0.8) evolved from the harmonic 
oscillator ground state. We contrast the non-resonant case, where 
rf = (pTi {(f> = {\ + V5)/2 is the golden ratio] for (a) j = 36 kicks 
and (b) j = 108 kicks, with the resonant case q^ = n for (c) j = 36 
and (d) j = 108. Although Qj(a,a') > 0.04 close to the origin in 
(b), our scale is fixed at this maximum to aid comparison between 
the images. 



The final expression obtained for M„ „ (29 1 is identical to 



that for M,'^], ( |24[ ) but with ^ consistently replaced by 2^. Fur- 
ther evolution will see such cycles repeated; the time evolution 
of the coefficients is entirely contained within the changing ar- 
guments of the Bessel functions: 



(30) 



where 



/ _ 



2(N-2)/2 

2(N-y)/2 ^ ojjjj 
2(N-V>/2 ^odd 



(31) 
(32) 



As shown in Fig. [T] in the case of quantum resonance the 
phases of the lattice state coefficients [Eq. ( 30 1] take an ex- 
tremely simple form, which repeats regularly across all of 
phase space. We emphasize that this holds for any initial co- 
herent state. 

We have chosen the case rf' - n for simplicity, however the 
above analysis holds for all 77^ - wn. For more general //, each 
application of Graf's theorem introduces a A:-dependent phase 
(see Appendix [B), preventing an end result as straightforward 
as (30 1. As a consequence, no such regularity of the phases, 
as shown in Fig. [l] is apparent for non-resonant values of rf 
beyond the first two kicks. 



4. Analytic solution at quantum resonance when q = 3,6 

In these cases the locations in phase space of the coherent 
states making up the lattice states ( [T4| are defined in terms of 
/ and ie^'^"^'^, which represent two non-orthogonal directions 
in the complex plane. This non-oithogonality is necessary to 
deiive the coefficient mapping (20i, but it adds complexity 



to the evolution of the coefficients for a coherent state initial 
condition. In contrast to the case q - A, where all coefficients 
M\i}„ are a direct function of m and n at quantum resonance, 
in the cases q - 1> and q - 6 the coefficients m]1]„ for y > 3 
can only be expressed as a summation over products of three 
Bessel functions, even at quantum resonance. 

Despite this added complexity, the cases <? = 3 and q - 6 
remain tractable at quantum resonance. There is a three-step 
cycle in the coefficients, analogous to the two-step cycle in 
the coefficients occurring when q = 4, with the coefficients 



M: 



IMj+i)] i^gjjjg equal to mI^JJ but with the arguments of all 

6 is 



Bessel functions doubled. A derivation for the case q 
given in Appendix [C] 



IV. MANIFESTATION OF QUANTUM RESONANCES 

A. Ballistic phase space expansion 

As for the atom-optical 5-kicked rotor (S) |46l or 5-kicked 
accelerator |50 5T\ systems, the primary observable effect of 
the quantum resonances in the KHO is an enhanced absorp- 
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FIG. 3: (color online). Mean energies yielded by KHO evolution 
when q = 4, K = -0.8, and the initial condition is the harmonic 
oscillator ground state, (a) Plots of the number of kicks required to 
reach a mean energy of 50% (lower line) and IOOHll) (upper line) as a 
function of rf. (b) Mean energy as a function of A' for rf = n, n/2, 
V3;r/2, and = (I + V5)/2 is the golden ratio]. 



tion of energy from the kicking potential. This is visible as 
an enhanced rate of expansion of the oscillator's state through 
phase space, as can be seen in Fig.|2] which shows Q functions 
(Husimi distributions) [481 Qj{a,a*) = \{iffj\a}f-/7: evolved 
from an initial harmonic oscillator ground state under resonant 
and non-resonant conditions. The evident symmetry between 
the position and momentum quadratures of the Q functions 
in Fig.|2]mean in situ measurement of the position density or 
time-of-flight determination of the momentum density would 
both provide equally effective signatures of quantum resonant 
evolution, and the Q function itself may also in principle be 
determined directly |29, 52 1. Note, however, that in any ex- 
periment the harmonic potential cannot be infinite in extent; 
after a large number of kicks the state of the oscillator will 
extend beyond the harmonic region. Consequently, the subse- 
quent time evolution will depart from the predictions of our, 
entirely harmonic, model. 

In the KHO, the expansion through phase space changes 
from diffusive (energy growing linearly in time) to ballistic 
(energy growing quadratically in time) at the quantum reso- 
nance values. In Fig.|3ja) we have plotted the number of kicks 
required for the KHO to achieve particular mean energy val- 
ues, as a function of jf when starting from an initial ground 
state. This reveals distinct minima in the number of kicks re- 
quired when jf = W7T, i.e., a quantum resonant value. Other, 
less prominent minima are associated with higher-order reso- 
nances, i.e., when jf is a rational multiple of n, with the rate of 
energy absorption falling rapidly as the rational denominator 
increases. Figure [Sj^b) shows explicitly how the mean energy 
grows quadratically as a function of the number of kicks when 
if - 71 ox njl (ballistic expansion), quite distinct to the two 
cases we have shown when if is an irrational multiple of n 
(diffusive expansion). 



spectrum of the KHO is a Hofstadter butterfly-type multifrac- 
tal ll38l l53l . as shown in Fig.|4] The theory of butterfly-type 
spectra in kicked quantum maps has been extensively studied 
||10l|4T]|42]|43]|44l|45]|54l: In summary, a pure point spectral 
component is associated with localized behaviour, whilst sin- 
gularly continuous and absolutely continuous spectral com- 
ponents are associated with diffusive and ballistic expansion 
respectively |43, 44|. For the KHO (when q e {3,4,6)) this 
means that the quasienergy spectrum of F switches between 
being singularly continuous at irrational multiples, and abso- 
lutely continuous at rational multiples of the primary quantum 
resonant values of rf-. 

Kicked atom-optical systems have been proposed as can- 
didates for direct observation of butterfly-type spectra |55|, 
and the atom-optical KHO considered here presents another 
system in which this might be realised. The KHO also of- 
fers the opportunity to directly explore the changes in the 
nature of the butterfly spectrum across the quantum reso- 
nance/an tiresonance transition (that is, as a function of the rel- 
ative alignment of the kicking and trapping potentials). These 
changes have been explored theoretically for the case q - A 



C. Experimental realization 

To observe, for example, the principal quantum resonance 
in the case where q - A requires a Lamb-Dicke parameter 
of T] = K y/H/lMo) = V^, and hence a suitable combina- 
tion of harmonic trapping angular frequency co and wavevec- 
tor K. This could be realized in a linear ion trap with rela- 
tively weak axial confinement, perhaps in conjunction with a 
standing wave formed by two counterpropagating laser beams 
overlapping at an acute angle. Typical experimental configu- 
rations involving laser cooled alkali atoms within a magnetic 
trap yield much larger values of rj; a configuration equivalent 
to that used by Duffy et al. would require a trap frequency 
of 4.8 kHz |24|. Such frequencies are nevertheless within 
reach of recent generations of magnetic chip traps |56|. Alter- 
natively, a potentially quite flexible set-up could be achieved 
by generating the kicking field with a CO2 laser (wavelength 
A - 10.6/im). Applying the kicking field directly along the 
axis of harmonic confinement, values of the Lamb-Dicke pa- 
rameter around the principal quantum resonances (e.g., from 
rf - njl up to rf- - 2n) could be reached using trap fre- 
quencies ranging from 8-36 Hz ('^^Cs) up to 40-200 Hz 
(^■'Na), and introducing the kicking field at an angle to the trap 
axis would allow these frequencies to be increased. Kicking 
strengths would be low because of the extreme detuning of the 
CO2 laser, but the quantum resonances will still be manifest 
at any value of k after a sufficient number of kicks. 



B. Hofstadter butterfly spectrum 

The quantum resonances are also apparent in the 
quasienergy spectra of F, i.e., the values ip such that e'f is 
an eigenvalue of F. Taken as a function of rf, the quasienergy 



V. CONCLUSIONS 

In this paper we have considered the characteristic quantum 
resonances occurring in an atom-optical 5-kicked harmonic 
oscillator (KHO), a system in which typically alkali metal 
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(a) (b) 




ff/n 

FIG. 4: (color online). Numerically determined quasienergy spectra of F, computed in a truncated Fock basis of 2000 states and plotted against 
jf. (a) and (b) show the overall structure (inset) and detail from the uppermost band in the cases q = A and q = 6, respectively. A Hofstadter 
butterfly-type structure f38"531 is seen in each individual band, the cell boundaries (points at which the spectrum becomes continuous across 
the band) of which indicate the quantum resonant values of rj^ = wn. One half of a cell for the case ^ = 4 is shown in detail in (c), better 
illustrating the division of the band into b continuous segments at rational fractions alb of the cell width. In each case k = -0.8, and eigenvalues 
are colored according to the overlap of the associated eigenstate with the harmonic oscillator ground state, from blue (~ 10"^^) to red (~ 10"'). 



atom(s) or an alkaline earth ion are trapped in an approxi- 
mately one-dimensional harmonic potential and periodically 
driven by pulses of an ofF-resonant laser standing wave. We 
have demonstrated that the atom-optical KHO displays quan- 
tum resonances dependent on the value of if, where rj is the 
Lamb-Dicke parameter, equal to the ratio of the width of the 
ground state of the harmonic trapping potential to the wave- 
length of the laser standing wave. 

When the oscillator is kicked a rational number of times rjq 
per oscillator period, with q equal to 3, 4 or 6, quantum reso- 
nances occur . At these values of q the eigenstates of the Flo- 
quet operator also possess a crystal phase space symmetry. 
We have shown that the quantum resonances occur at natural 
number multiples of a primary resonance value: 77^ = tt in the 
case q - A, and if - Inj V3 in the cases q - 3 and ^ = 6. At 
these values of if' all q kicks within commute. Further, sec- 
ondary, quantum resonances occur for if equal to non-integer 
rational multiples of the primary values; at these values of rf 
there is a partial commutativity between the kicks constituting 
F'i. We have demonstrated that the quantum resonances can 
be interpreted in terms of the relation between and its own 
eigenstates; specifically, the quantum resonances occur when 
F'' can be decomposed in terms of a set of displacement oper- 



ators which is equivalent to the set of crystal symmetries of the 
eigenstates. We have also derived an expression for the time 
evolution of a coherent state in the cases <7 = 3, 4, 6 by recast- 
ing the dynamics in terms of a mapping on a two-dimensional 
lattice of coefficients. At quantum resonance values of if the 
evolution of these coefficients is greatly simplified; we have 
shown that in this case there is a direct expression for all the 
coefficients at a given step. 

Finally, we have considered the manifestation of quantum 
resonances in an atom-optical KHO, in particular the ballistic 
increase in the mean energy associated with quantum reso- 
nances, and the quasienergy spectrum. To this end we have 
also discussed experimentally reasonable parameter regimes 
in which to observe quantum resonant dynamics in the atom- 
optical 5-kicked harmonic oscillator. 
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FIG. 5: (color online). Pictorial representation of the values of 
sm{2nj/q) (the vertical position of points) in the cases q = 3, 4, 
5, and 6. 



(d) z„ = V3/2 for ^ = 6 (tj^ = win/ ^/3 at quantum reso- 
nance). Figure [Sj^d) also illustrates the absence of quantum 
resonances for q > 6; the addition of another distinct point on 
the circumference of the circle shown necessitates at least a 
second value of z„. 



for support. 



APPENDIX B: GRAF'S ADDITION THEOREM 



APPENDIX A: QUANTUM RESONANT VALUES OF r/^ 

Because we are free to choose the sign of w on the right 
side of (|9]l we can, without loss of generality, simplify that 
commutation condition to 

77^|sin(2;r;/^)| = w;r; ; e [0, . . . , ^ - 1], w = 0, 1, 2, . . . . 

(Al) 

Quantum resonances occur when the condition ( |Al| l is satis- 
fied for all possible j. 

Ignoring the trivial cases q - 1 and q - 2, where 
smilnjiq) - for all possible j, | sin (27rj/^) | takes either 
iq - 2)/2 values other than zero (even q), or {q - l)/2 values 
other than zero (odd q). If we denote these nonzero values 
by Zn [where n ranges from Xo (q — 2)/2 (even q), or from 
to (^ - l)/2 (odd q)], the condition for quantum resonance 
becomes 



(A2) 



Because the z„ are the moduli of sines equally distributed 
in angle, they cannot share a common factor other than unity. 



Consequently, condition ( A2 1 can be satisfied for all z„ only 
if all of the z„ are equal. The cases of interest are enumerated 
pictorially in Figure pi where we see that: (a) z„ — V3/2 for 
q - 3 {rf — win I at quantum resonance); (b) z„ = 1 for 
q = A {if - wn at quantum resonance); (c) there is more than 
one value of z„ for ^ = 5 (there are no quantum resonances); 



Graf's addition theorem for two Bessel functions of integer 
order and equal, real argument ^ is ISTl 



(Bl) 



/c=-oo 



where (' and x are related to ^ and a by 



^' = ^ V2[l - cos(a')], (B2) 
ai-cos(a)] =^'cosOr), (B3) 
^sin(a) = ^'sinCr). (B4) 

These relations are shown graphically for positive f in Fig. |6] 
for < Q' < ;7r. In the quantum resonant case we consider in 
section [niC3| a - k. Consequently = 0, and we obtain 



2 Jn^uiOJkiOe""' = Jn{20. 



(B5) 



APPENDIX C: EVOLUTION OF AN INITIAL COHERENT 
STATE WHEN g = 6 

In the case q - 6, and at quantum resonance rj^ - 2n/ V3 



the mapping ( 20 1 becomes 

oo 

Mlc:' = 2 i'm)Mi^„_,^_j^-^'"-. (CI) 
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As in the case = 4, we begin with the single coherent state 
initial condition 



Replacing the indices in the correct fashion we get 



[0] 



m+n-k.-tn 



and thus, substituting this into the mapping (CI 



(C2) 
(C3) 
(C4) 



Again, as in the case ^ = 4, we observe the more general 
pattern of evolution by repeating this process u ntil we obtain 



the coefficients M|,f |,. We begin by rearranging (|C4|i to get 

[1] 



-i{m+n-k)in-k)n 



(C5) 



Substituting this into the mapping (CI i produces 

M]^„ = r^"Jr„UOi"'Jm{Oe-''"^"'^"^\ (C6) 
where we have used the identity /"'V_„(0 - i"Jn(0- Rear- 



ranging ( C6 1 we find that 



and hence 

<i = X i'Jk{0i"-''Jn-k{Oi'"^"-'Jm.n-k{Oe-' 



(C7) 



-i{mn+n- +k~)n 



k=-c 



(C8) 

Graf's addition theorem (Appendix |B]i cannot be directly 
applied to ( |C8| l, as there are three Bessel functions with k- 
dependent indices. However, we accept this summation as the 
final form of M^j, and continue with the process, using k' as 



the new index of summation. Rearranging ( C8 1 we have 

oo 



-k-k' 

k=-oo (C9) 

i[(m+n-k' )m-m^-k-]7: 



XJn-k-k'iOe'^ 



which, when substituted into the mapping (CI i, produces 

oo oo 



k=-oo k'=--oo 



J(mn-k-)n 



Z:k J / >-\:m+k j / i-\;n-k i{nm-k^)n 
I JkiO' Jm+k(i)l 

OO 
OO 

^ JAOJn-k^k'iOe"'''', 



k=- 
X 



(CIO) 



k'=-co 



= 2 i''m)i"'^''J„,^k(0i"'''Jn-k(20e'^""'-'"^\ 

k=-co 

where we have used Graf's theorem to ehminate the summa- 
tion over k'. Rearranging this we obtain 



m+n-k',-m 



■m+n-k'+k j / 
* ^m+n-k'+k\^) 



(Cll) 



-i[m(m+n-k' )+k^]7T 



and hence 



oo oo 



M'i = Z Z i''MOi'jk(Oi'"^"-''^'j,n.n-k'.k(0 

k--oo k'--oo 



= 2 i''m)i'"^''J„,u(20i"'^"^''e-'^ 



\m{m+n)+k^]ji 



k=- 
X 



^ Jk'iOJ m+n+k' +k (0^''' " 



(C12) 



k'=-o 



^:=-oo 



Repeating the process one final time we have 



M^^L-k'^-^n = Z i'M0i"''"-''''jn..n-k'.ki20 



-i(mn+rp-+k' +k^-mk')K 



and hence 



oo oo 



= Z Z i''M0i'jki0i"'^"-''^'jn,.n-k'.ki20 



A:--oo k'--c> 



X I Jn-k'+km)e 



If we define k - k' -\- A, then < |C13[ ) takes the form 



<!, = Z ''""""^-/-«+«+a(20/""^^«+a(2^) 



(C13) 



A=-oo 



(C14) 



^'=-00 



By applying Graf's theorem to the summation over k' in 



(C14i, and then applying the transformation A — > -A we ob- 
tain 



Mini = Z '^M20i"'^J.,-A(20 



(C15) 



X /'«+"-^7,„+„_A(2^)e 



that is, an analogue of M„ „ with f replaced by 2^ and k re- 
placed by A. As in the case q - 4, the evolution from this 
point on will continue cyclically. 



-i{mn+n +A-)7i 
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